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ABSTRACT 

A star formation efficiency per free fall time that evolves over the life time of giant molecular clouds 
(GMCs) may have important implications for models of supersonic turbulence in molecular clouds or 
for the relation between star formation rate and H 2 surface density. We discuss observational data 
that could be interpreted as evidence of such a time variability. In particular, we investigate a recent 
claim based on measurements of H 2 and stellar masses in individual GMCs. We show that this claim 
depends crucially on the assumption that H 2 masses do not evolve over the life times of GMCs. We 
exemplify our findings with a simple toy model that uses a constant star formation efficiency and, yet, 
is able to explain the observational data. 

Subject headings: galaxies: evolution — stars: formation 



1. INTRODUCTION 

The life times of giant molecular clouds (GMCs) have 
been at t he center of a major deba t e for at least the last 
40 yea rs (iGoldreich fe Kwarj||1974t IZuckerman fe Evansl 



[19741: iSolomon et al.1 Il979t lElmegreeri I2000L 120071) 

GMCs that live for many free-fall times need a 
mechanism that prevents them from gravitational col- 
lapse. Over the last couple of years the consensus 
is growing that the life times of GMCs are likely 
a few free-fall times, or even l ess (Q jlmegrccn 120001 : 
iBallesteros-Paredes fe Hartmannl l2007t iMurrlrvl 120101 ) 
and the focus has shifted towards the challenge of ex- 
plaining the low star formation efficiencies in GMCs. The 
star formation efficiency per free fall time eg is defined 
as the ratio of free fall time tg to gas depletion time 
Mh 2 /M*. In other words: 



iff 



(1) 



i.e. the instantaneous star formation rate (SFR) is pro- 
portional to the available amount of molecular hydrogen 
(H 2 ) via the proportionality fac tor eg/tg. The observed 
value eg - 0.01 - 0.02 (e.g. iKrumholz fc Tanl I2007H 
means that only 1-2% of the mass of a GMC is converted 
into stars over a free-fall time. If s tar formation is sup- 
porte d by supersonic turbulence ([Krumholz fc McKed 
120051 ) eg is expected to be only very weakly dependent on 
the Mach number of the turbulent flow and thus approxi- 
matively constant, but this may be an oversimplification, 
see e .g. ([Vazquez- Semadeni et al.ll2005l : iLi fc Nakamural 
I2006T) . On the other hand, if GMCs have life times of 
the order of a free fall time they do not need be sup- 
ported by turbulence. The star formation efficiency in 
such clouds may increase as the Mach num ber in the flow 
decre ases and the cloud collapses, see e.g. iBonnell et al.l 
120101 We note that a time varying eg should introduce 
additional scatter in the relation between star formation 
rate and H 2 surface density on small (< 100 pc) scales. 
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This scatter should propa gate up to ~ kpc scales, see 
e.g. ([Feldmann et alJl2010D . and hence would contribute 
to the scatter in the Kennicutt-Schmidt relation. This, 
at least in principle, could be used to test observationally 
the time dependence of eg. 

In section [5] we will discuss two common misconcep- 
tions that could give rise to the impression that eg varies 
over the life times of GMCs even if it is a constant. In 
section [3] we present and analyze a toy model in order to 
exemplify and quantify our statements. 

2. DO OBSERVATIONS CONFIRM A TIME- VARYING STAR 
FORMATION EFFICIENCY? 

In equation ([I]) we define the star formation efficiency 
per free fall time eg. Another commonly used efficiency 
is the star formation efficiency of the GMC cgmCj i-e. 
the fraction of H 2 mass of the cloud that is converted 
into stars over the life time of the cloud. In a picture 
where GMCs start with an initial reservoir of H 2 , which 
is used in the subsequent star formation process, the final 
stellar mass (final) is divided by the initial H 2 mass of 
the cloud. If the cloud accretes a substantial amount of 
H 2 over its life time, the definition has to be generalized. 
We will use: 

M* (final) 



£gmc 



max(Affj 2 



where max(Mn 2 ) is the maximal H 2 mass of the GMC. 
By definition cqmc is a non evolving quantity and it can 
be estimated, e.g., by comparing luminosity distribution 
of OB associations in the Milky Way with the mass s pec- 
trum of molecular clouds ([Williams fc McKeejl"l997D . It 
cannot be directly measured on a cloud-to-cloud basis, 
because M* and Mh 2 must be known at two different 
times . Instead such observations, see e.g. (jMvers et al.1 
1986), estimate the following quantity 



r]GMc(t) = 



M„(f) 



( M«(t) 

M na (t)+M.(t) ~ M H2 (t)' 



(3) 



The latter, approximate equality is due to the fact that 
for most observed GMCs is smaller than Mh 2 . Obvi- 
ously, ?7gmc (t) increases over the life time of a cloud and 
should not be confused with either cqmc or eg. From 
([T]) we can estimate M*(final) = £eg/tg max(MH 2 )tfmai> 
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hence ecMC = £ e a^final/%) where £ is a constant fudge 
factor of order unity that depends on the actual time 
evolution of the SFRs and Afn 2 (and es if it is time de- 
pendent). We will estimate £ for a simple toy model in 
section [3l Combining this result with equation ([2]) and 
((3|) we obtain: 



VGMc(t) « £e ff 



^final 

~ts~ 



M* (final) 



max(Mnj ) 



M U2 (t) 



(4) 



There are several ways of creating large values of ?7gmc 
and they correspond to the various terms in this equa- 
tion. First, eg could be time dependent. For instance, it 
could smoothly increase as the cloud collapse advances 
or, alternatively, vary stochastically about some aver- 
age value. A second possibility is that some clouds may 
live for many free fall times, i.e. tfinaiAff is large in a 
subset of GMCs. The third factor in the third bracket 
in equation explain why ?7gmc can also be smaller 
than £gmc- Finally, ijqmc can be boosted if the ob- 
served H2 mass is significantly less than max(A/H 2 ), i.e. 
if GMCs lose (in one way or another) a large fraction of 
their molecular hydrogen over their life time. The lat- 
ter scenario predicts that t]gmc (t) should roughly scale 

oc Mt7 over the life time of individual GMCs. An ob- 
? 

servati onal sample o f an ensemble of GMCs shows this 
trend (Murray l2010f ). However, this trend can also be 
produced by a selection effect based on stellar mass, e.g. 
selecting GMCs with M* > M*ji m i t excludes values of 
jjgmc that are small e r tha n M* j ii m i t /MH 2 , see equation 
(J2J). In fact, iMurravl (|2010t ) is selecting clouds based on 
ionizing luminosities, which roughly corresponds to se- 
lecting clouds based the stellar mass formed within the 
last 4 Myr. Such a selection effect explains why a dif- 
ferent s tudy of ~ 10 5 M(7) GMCs find much lower effi- 
ciencies lLada et "all (|2010D . The existence of the selec- 
tion effect is not an argument against or in favor of an 
evolving eg, rather it shows th at the GMCs with large 
values of ?7gmc in the sample of lMurravl (|2010t ) are likely 
a heavily biased subset. In the case that eg is, in fact, 
a non evolving quantity and the measured large values 
of t?gmc are driven by changing molecular gas masses, 
we can make a rather generic prediction. The similar- 
ity of the scaling with GMC mass 4 (oc M£ ) of ^gmc, 
on the one hand, and the lower boundary of the region 
excluded by the discussed selection effect, on the other 
hand, implies that the observed GMCs with large values 
of ijgmc should have rather similar maximal H2 masses 
max(Mn 2 ). The toy model that we discuss in section [3] 
predicts max(M H2 ) ~ 10 6 - 1O 7 M . We note that this 
scenario explains rather naturally the absence of massive 
(> 10 6 M Q ) GMCs with high values of ijgmc- 

A different issue can arise if one compares star forma- 
tion rates and H2 masses in order to estimate eg/iff via 
equation (QJ. For example, let us assume that we mea- 
sure SFRs and H2 masses within small (< 100 pc) aper- 
tures around peaks of CO emission (tracing the H 2 mass) 
and peaks of Ha emission (tracing star formation rates) , 

4 A linear regr ession of ijgmc vs - -^h 2 for the data presented 
in Murray (2010) gives a slope of —0.59 ± 0.19. This is consistent 
with the prediction of our toy model (slope ~ — 0.75, see section 
[3j that takes into account that, in fact, not the total stellar mass 
has been measured, but only the stellar mass formed within the 
last ~ 4 Myr. 



see e.g. iSchruba et all (f2010h . If we observe that CO 
peaks have lower SFRs at given H2 mass compared with 
peaks of Ha emission, does this imply a time-varying 
effAff? The answer to that question depends on the way 
the SFRs are measured. SFRs that are derived from 
Ha emission are effectively averaged over the past 5-10 
Myr, which might well be a significant fraction of the life 
time of the molecular cloud. For SFRs that are based on 
Ha+24fj,m emission this averaging time span would be 
even longer. The star formation efficiencies per free-fall 
time that are estimated from such a time averaged SFR 
will be small initially (no stars have been formed over 
most of the time averaging interval simply because the 
GMC has only formed recently). The measured SFRs 
will increase until the age of the GMC is similar to the 
averaging time span. In addition, the H2 mass of the 
cloud might evolve (possibly decrease) leading to an ad- 
ditional increase in the apparent value of eg/tg with time. 
If the following three conditions are satisfied, a difference 
in the measured SFR per measured H2 mass can provide 
strong evidence for a time-varying star formation effi- 
ciency per free fall time. First, the averaging times of 
the SFRs need to be small compared to ages of the ob- 
served clouds. Second, the observable H2 reservoirs need 
to be close to max(Afn 2 ), and, finally, the free fall times 
of the clouds need to be known. A recent study that 
meas ures SFRs with re asonably short averaging times (2 
Myr, lLada et"aLll2010D estimates star formation efficien- 
cies per free fall time of the order of 2% for most clouds 
in the sample, with the scatter mostly driven by the mass 
of molecular gas of relatively low density (n < 10 4 cm~ 3 ) 
that does not participate in the star formation. 

3. TOY MODEL 

We will now discuss a toy model in order to both ex- 
emplify the points made in section [5J but also to provide 
a framework in which we can make some quantitative 
predictions. We should stress that the statements made 
in the previous section are completely generic and do not 
depend on the specific assumption that go into the model 
that we are going to present. Our model is almost insult- 
ingly simple, and, given that, our aim is not to reproduce 
the full complexity in the evolution of GMCs or even, to 
be consistent with any available observation. On the 
other hand the model offers a pragmatic approach to the 
mass evolution of GMCs and may be easily generalized 
to facilitate more complex scenarios. 

The ansatz of the model is to supplement equation (JTJ) 
with an equivalent equation that describes the evolution 
of the H2 mass: 



Mn 2 = -?M Hj -aM,+7 
iff 



(5) 



The extra term aM, is motivated by assuming that stel- 
lar feedback is limiting the life time of molecular clouds, 
e.g. via photo-ionization, thermal pressure or radiation 
pressure (IWilliams fc McKeelll997t iMurrav et~aL1 l2010t 
lLopez et alJl2010T ). This feedback should therefore cou- 
ple to the formed stellar mass via some efficiency factor 
a that sets the time scale for the destruction/removal 
of H 2 from the cloud 5 . The term 7 is the net "accre- 

5 Depending on the type of feedback M* should refer to the total 
stellar mass times a weight parameter that takes into account that 
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tion" rate of H2, which includes all processes that cre- 
ate and destroy H2 and are not directly coupled to ei- 
ther M* or Mh 2 . Both a and 7 could in principle be 
time dependent. For simplicity we assume that they are 
cons tant. Our model is minimalistic (compared with , 
e.g..lMatznerl[2r)0llTan et alJl20M fluff fc Stahledl^fM 



Krumholz ct al. 2006), but it has the advantage that we 



can parametrize our ignorance of the relevant physical 
processes that destroy and disperse the cloud into the 
parameters a and 7. Together with appropriate initial 
conditions equations ([T]) and ([5]) fully determine the evo- 
lution of the masses of molecular hydrogen and the stellar 
component in a CMC. 

We will also make the simplifying assumption that the 
free-fall time does not evolve strongly over the history 
of the GMC, i.e. both the star formation efficiency per 
free fall time and the star formation time scale are now 
fixed. This assumption is not crucial for the model, but 
we will use it for the following reasons. First, there is no 
clear systematic trend of free fall time with mass over the 
ra nge of GMCs t hat we are comparing to, see e.g. Table 2 
of lMurravl(l20Tol) . Second, assuming a non-evolving free- 
fall time allows for a convenient analytical solution of the 
problem. Third, we find that even with this assumption 
our model describes the observed data reasonably well. 
We stress that our main aim is to show that a simple 
model can produce an observational signal that could 
be misinterpreted as evidence for evolution of the star 
formation efficiencies. We do not try to model the precise 
properties of the ensemble of GMCs in the Galaxy. 

With tg fixed (and, of course, we assume that the star 
formation efficiency per free fall time is a constant, too) 
we can insert ((T|) into ([5]) and obtain a linear 2nd or- 
der differential equation for Mjj 2 , i.e. the equation of a 
damped harmonic oscillator. 

Solving the differential equation we obtain 



Mh 2 (t) = Ae~ tb/2 cos(ujt + 
M H 



M„(i): 



(wtan(wi + 0) - 6/2) 



7 



(6) 
(7) 



a a 

where b — eg/tg is the inverse of the star formation 
timescale, and ui = \J ab — b 2 /4 is the "oscillation" pe- 
riod. 

Phase </> and amplitude A depend on the initial condi- 
tions. In the following we restrict ourselves to two special 
cases of the general model ©, ((?]). 

• No accretion scenario: Assumes 7 = 0, Mfj Q (i — 
0) = M > 0, and M*(i = 0) = 0. It follows 
(f> = atan(6/(2w)), and A — M / cos(4>). 

• Pure accretion scenario: Assumes that all H 2 is 
"accreted", i.e. Mjj 2 (t = 0) = 0, M*(f = 0) = 
and 7 > 0. In this case phase and amplitude are 
given by = — 7r/2, A = j/lj. 

We adopt the parameters eg = 0.02 and tg = 6 Myr, 
which are consist ent with observations of eff over a range 
of density scales (|Krumholz fc Tanll2007l ). and with the 
free fall times 6.1^40 Myr measured in the sample of 

feedback is provided by stars which have a limited life time. For 
simplicity we will assume that M* is the total amount of stellar 
mass formed within the cloud. 



iMurravl (|2010|) , respectively. We note that only the ratio 
eg/tg = 0.0033 Myr -1 enters our model. The a param- 
eter is chosen such that the life time of the cloud, i.e. 
the time tg na i at which MH 2 (tfmai) = 0, is ~ 20 Myr 
([Williams fc McKed [19971 ) . Hence, we use a = 2Myr~ 1 
in the no accretion scenario and a = 8 Myr - 1 in the pure 
accretion scenario, respectively. 

Assuming eg/tg <C a the life time of a GMC is given 
by 



^final 



, and tg 



inal 



y/ae s /t s 



The left (right) expression refers to the no accretion (pure 
accretion) scenario. We note that in both considered sce- 
narios the life time does not depend on the initial cloud 
mass or the accretion rate, respectively. The evolution 
of Mh 2 , M* and M* :< 4M yr normalized to max(Mn 2 ) is 
shown in Fig. [TJ 

Assuming eg/tg -C a we can easily estimate the total 
stellar mass that is formed during the life time of the 
cloud from equation (JJJ . In the no accretion scenario we 
obtain 



M*(t fina i) « M 



eff Aff 



f 



ifinal , 
Ztff 



Mo 



effAff 



a 



while the pure accretion scenario predicts 
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M*(tg na i) « 

a 



ifinal , 

Ms 
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In the pure accretion scenario a GMC attains its maxi- 
mum mass at t w tfmal/2- The H2 mass is then approx- 
imatively 7 / 1 yjaegjtg. Combining these results we see 
that the star formation efficiency of a GMC is 



cgmc 



eff Afl 



, and e G MC ~ 2 



effAfi 



Again, the left (right) expression refers to the no ac- 
cretion (pure accretion) scenario. Written in terms of 
the life time of the GMC both expression are identical, 
namely e G McAfinai ~ ( 2 /tt) eg/tg, i.e. £ = 2/n (section 

In Fig. [2] we show the predictions for 7/gmc and 
rjg of the two scenarios of our model, toget her with 
eff an d £gmc, and the observational data from IMurravl 
(2010). To be consistent with observations only the stel- 
lar mass that formed within the last 4 Myr is included 
in the definition jjgmc, see ((3]), and rjg is estimated as 
Vs = Vgmc tg/4 Myr. Our model reproduces the trends 
of tjgmc and rjg with GMC mass, suggesting that these 
are maybe not solely due to selection effects. In both 
scenarios j/gmc and rjg roughly scale as M^ 0,75 over the 
mass range 10 4 M Q - 10 6 M Q . 

With the chosen parameters our model predicts that 
Vgmc is only significantly larger than 6gmc for the last 
~ 1 Myr in the life of a GMC, this includes most of the 
G MCs with mas ses less than ~ 10 6 M Q in the sample 
of IMurravl (j2010T j. We note that the precise time does 
depend on the assumed life time of the cloud. Clouds 
with shorter life times spend more time in a state in 
which jjgmc > e GM c- 

Our model exemplifies that it is difficult to prove 
the existence of a time- varying star formation efficiency 
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Fig. 1. — The evolution of the masses of the GMC components (normalized to the maximum H2 mass of the GMC) according to the 
two scenarios: no accretion (left) and pure accretion (right). We assume eg/tg = 0.0033 Myr -1 and ct = 2 Myr - 1 (a = 8Myr -1 ) in the 
no accretion (pure accretion) scenario. The different lines correspond to the H2 mass (solid blue line), total stellar mass (dashed red line), 
stellar mass formed within 4 Myr (dot-dashed green line). 
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Fig. 2. — Estimators of the star formation efficiencies as function of mol ecula r mass of the GMC . The quantity ijgmc (Vff) is shown in 
the left (right) panel. The observational data presented by Murray (2010) and Lada ct al. (2010) are indicated with empty squares and 
triangles, respectively. Overplotted are solid (dot-dashed) lines that refer to the predictions of the toy model in the no accretion (pure 
accretion) scenario. The 6 solid and dot-dashed blue lines at the top use en/iff = 0.0033 Myr -1 , a = 2Myr _1 (ct = 8 Myr - 1 ) in the no 
accretion (pure accretion) scenario and should be compared to the squares. Note that for consistency only the stellar mass formed within 
the last 4 Myr is considered in the computation of jjgmc an d ??ff. Red lines (near the bottom left of each panel) use free fall times that 
are a factor 2 smaller (since these are intrinsically smaller clouds), a that are a factor 2 larger (to keep the same cqmc) an d the stellar 
mass that formed within the past 2 Myr (the masses are derived from counting young stellar objects). These lines should be compared to 
the triangles. Filled circles and filled stars indicate when the age of the modeled GMC is half its total life time or when the cloud is 1 Myr 
away from the end of its life, respectively. The diagonal dashed line indicates a slope of -0.75, which is approximatively the slope predicted 
by our toy model. The observed data is consistent with this slope. A linear regression of ??gmc an d TJs as function of GMC mass using all 
clouds with masses > 10 4 A/q returns slopes of —0.59 ± 0.19 and —0.49 ± 0.32, respectively, at the 95% confidence limit. Each panel also 
contains a horizontal line that denotes the value of the star formation efficiencies eg = 0.02 and cqmc = 0.04, respectively. 
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based on observational quantities such as ?7gmc or Vfi- We thank A. Kravtsov, M. Krumholz and N. Murray 

This is not to say that such a time-dependence does not for helpful comments, 
exist, we merely conclude that current observational ev- 
idence for its existence is insufficient. 
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